High Energy Neutrinos from Novae in Symbiotic Binaries: The Case of V407 Cygni 
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Detection of high-energy (> 100 MeV) 7 rays by the Fermi Large Area Telescope (LAT) from 
a nova in the symbiotic binary system V407 Cygni has opened possibility of high-energy neutrino 
detection from this type of sources. Thermonuclear explosion on the white dwarf surface sets off 
a nova shell in motion that expands and slows down in a dense surrounding medium provided by 
the red giant companion. Particles are accelerated in the shocks of the shell, and interact with 
surrounding medium to produce observed 7 rays. We show that proton-proton interaction, which is 
most likely responsible for producing 7 rays via neutral pion decay, produces > 0.1 GeV neutrinos 
that can be detected by the current and future experiments at > 10 GeV. 

PACS numbers: 95.85.Ry, 98.70.Sa, 14.60.Pq 



I. INTRODUCTION 



High-energy (> 1 GeV) neutrinos are produced dom- 
inantly via decays of charged pions and kaons created 
by proton-proton (pp) and proton-photon (jry) interac- 
tions. Powerful astrophysical sources of 7-rays, long du- 
ration (> 2 s) gamma ray bursts (GRBs); short duration 
(< 2 s) GRBs and active galactic nuclei (AGNs), as well 
as sources such as core-collapse supernovae, supernova 
remnants and microquasars have been proposed as the 
sources of high-energy i/'s (see e.g. [l]). Protons and ions 
accelerated by a Fermi mechanism in the shocks of these 
sources interact with ambient particles and/or soft pho- 
tons to produce 7 rays, via neutral pion decay, and i/'s. 
Modeling 7 ray emission from an astrophysical source by 
a 7T° model thus inevitably predicts high-energy v flux 
from the same source (see e.g. |2j). Detection of these 
i/'s can provide a conclusive proof of the ir° model and 
discriminate against a leptonic model for observed 7 rays. 

Detection of > 100 MeV 7 rays by the Fermi Large 
Area Telescope (LAT) from the Nova 2010 on March 10 
in the symbiotic binary V407 Cygni is the first from 
any nova |3]. The binary system consists of a Mira- 
type pulsating star, a red giant (RG) which may be 
swollen to a radius ~ 500i?p), and a white dwarf (WD) 
with mass Mwd > 1 Mq [4j . The WD accretes mate- 
rial from the stellar wind of the RG and forms an en- 
velope on its surface. As the mass of the envelope in- 
creases, an increasing WD surface temperature ignites 
a thermonuclear runaway (see e.g. Q). A shell of 
mass M e j can be ejected from the WD surface when the 
proper pressure at the core-envelope interface exceeds ~ 
2 • 10 20 (Af cj /10- 6 M Q )(A% D /M Q )(i?wD/10 8 cm)- 4 dyn 
cm~ 2 . Note, however, that the WD mass and radius are 
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related roughly as i?wD oc M wr j . Thus a smaller enve- 
lope mass and corresponding M j is required to produce 
a thermonuclear explosion on the WD surface as M\yD 
increases. 

An optical nova is detected when the initially hot shell 
expands to a large radius from the WD, thus cooled to 
a temperature < 10 4 K, and when the opacity dropped 
due to recombination of ionized hydrogen gas Non- 
thermal 7 rays detected by LAT for 15 days following 
the optical detection requires particle acceleration and 
their interactions in the nova shell and in the surrounding 
medium of the RG. Particle acceleration in the shock 
of the nova shell has been perceived for the 2006 nova 
outburst in another symbiotic binary RS Ophiuchi ja], 
but LAT detection of > 100 MeV 7 rays from V407 Cygni 
prove existence of these high-energy particles in the nova 
shock and their interactions. 

In this brief report we calculate the expected v flux 
from the nova 2010 in the symbiotic binary V407 Cygni 
following the ir° decay model of 7 ray emission from pp 
interactions reported in Ref. [3]. We also estimate ex- 
pected v events in current and future detectors and dis- 
cuss detection prospects of these type of novae. 



II. BINARY SYSTEM AND PARTICLE 
ACCELERATION IN NOVA SHELL 

Modeling of 7-ray light curve from the nova 2010 in 
V407 Cygni system suggests a binary separation of a s» 
10 14 ai4 cm, a mass-loss rate of M w s» 3 • 10~ 7 A+M Q yr^ 1 
from the RG that blows a wind with a measured velocity 
of v w « 10d Wi 6 km s _1 and a shell kinetic energy of ~ 
10 44 erg with a measured velocity t> j = 3200±345 km s _1 
[3J. Close proximity between the WD and RG allows 
the parts of the nova shell that expand toward the RG 
to decelerate and enter the Sedov- Taylor phase rapidly, 
resulting in efficient conversion of shell kinetic energy into 
shock-accelerated particle energy. 
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For a distance R and polar angle 9 from the WD cen- 
ter towards the RG, the density of particles in the RG 
wind can be calculated as n(R,6) — (M w /v w m)(R 2 + 
a 2 — 2aRcos9)~ 1 using an iverse square-law density pro- 
file of the RG wind. The magnetic field in the forward 
shock that propagates into the surrounding medium with 
a velocity u s h ~ u c j ~ 3.2 • 10 8 w s h,s.5 cm s -1 , before sig- 
nificant deceleration, can be calculated from the ther- 
mal energy density in the RG wind with temperature 
T w a 700T W o 8 K 1 as 



B(R,i 



0.04 



8kT w M w /v w m 
R 2 + a 2 -2aRcos9 
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G; 9 = 0, 



(1) 



in t = td day time scale. Here in — 10 _ g is the mean 
particle mass; is a magnetic field amplification fac- 
tor, which may arise from shock modification, and we 
assumed a factor 4 increase in particle density expected 
to arise in strong shocks. The maximum proton energy, 
accelerated by a Fermi mechanism, in this magnetic field 
can be estimated by assuming that the acceleration time 
scale is equal to the Larmor time scale as 



E, 



eBtvi 
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GeV; 6» = 0. 



(2) 



Here we assumed that tp — 20<y9i.3 is the number of e- 
folding required for a thermal proton with ~ 1 keV ki- 
netic energy to reach E praax . For 9 = n, i.e. in the part 



of the shell opposite to the RG, E P:Uiax 
t = t d day. 

Evolution of the shock velocity and radius in the de- 
celeration phase at later time (> 3 day for 9 ~ 0), af- 
ter the shell accumulates swept up RG wind material of 
mass equal to its own in a given solid angle, depends on 
the detail calculation for different 9 and is beyond the 
scope of this paper. Here we simply assume that pro- 
tons can be accelerated to an energy up to the value in 
Eq. ([2]) throughout the 15 day 7 ray emission episode, as 
required by the 7r° model. 



III. NEUTRINO FLUXES ON EARTH 

The observed 7 ray flux, averaged over 15 day outburst, 
from the nova in V407 Cygni at a distance D is fitted with 



a 7r — !• 27 decay model 



as 



^(E^^fjLtM Q(E P , E 1 )Np(E p )dE p , (3) 

Here Q{E p ,E 1 ) is the 7-ray production rate per unit 
density of H atoms Q, £m = 1-84 is the nuclear enhance- 
ment factor to take into account the contribution of other 
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FIG. 1: Neutrino source fluxes (thick curves) from Nova 2010 
in the V407 Cygni symbiotic binary calculated, using Eq. @, 
from the best-fit tt° decay model parameters (thin solid curve) 
of 7 ray emission. The data points are taken from Ref. Q| 



atoms Q, N p (E p ) — N p fiE~ s exp[(m p c 2 — E p )/E c , p ] is 
the shock-accelerated proton spectrum and (uh) is the 
average number density of target particles. The best-fit 
parameter values are s = 2.15^q 28: E C . P = 32lg 5 GeV 
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with a total energy in protons ss 6.9 
10 8 cm _3 ) _1 (L)/2.7 kpc) 2 erg in steadv 
that a delta-function approximation [9J 
(c(n H )/4nD 2 K w ) f~ dE^l - m 2 c 4 )" 



10 42 ((n ff )/4. 



state 



Note 



as $ 7 (-E 7 ) 
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90 GeV at E v f K 7! )N p {m p (r + E^/K^) also produces reasonable fit 
to the observed 7-ray flux with the same spectral pa- 
rameters s and E C:P . Here E n ,tb — E 1 + m 2 c A /AE 1 is 
the threshold pion energy and w 0.17 is the mean 
fraction of the proton kinetic energy converting to 7r°. 

Neutrino fluxes from -n^ decays, which are created in 
the pp interactions along with ir° , can be estimated from 
the observed 7-ray source flux as (Tol - [l2| 



$*J«(E) = = 0.50$ 7 (£); 

\E) = 0.30$ 7 (£); <5> S ^{E) = 0.22$ 7 (£), (4) 



for s — 2.0. For softer indices, s > 2, of the proton spec- 
trum the v fluxes decrease compared to the 7-ray flux, 
following the hadronic cascade theory [To[ . Contributions 
from kaons to the neutrino fluxes are at the level of 10% 
or less and we ignore those. 

Figure [1] shows the v fluxes (thick curves) calculated 
using Eq. (U), and the ir° decay model (thin solid curve) 
for the 7 ray spectrum from Ref. Q . The data points and 
2-er upper limits are also taken from Ref. Q. Note that 
the source fluxes will be modified as they are measured 
by any Earth-based detector due to flavor oscillation in 
vacuum and in matter (MSW effect) inside the Earth, 
which is particularly important below ~ 50 GeV. 
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High-energy v's, created with definite flavors (y a 
with a — e, /i, t), arrive as coherent mass eigen- 
states (vi with i = 1, 2, 3) from V407 Cygni on the 
Earth's surface. This is because the separation be- 
tween the wave packets of different mass eigenstates 
d s = LAm 2 /2El fa 2.4 • 10- 4 (i/2.7kpc)(Am§ 1 /2.4 • 
10 -3 eV 2 )(E 1/ /10 GeV)~ 2 cm, arising from propagation, 
is much smaller than the typical size of the v wave pack- 
ets from tt^~ and p^ decays. However, the v detec- 
tors do not have necessary energy resolution, AE/E ps 
10~ 13 ( J L/2.7kpc)- 1 (^/10GcV), to measure the oscilla- 
tory pattern arising from coherent v mass eigenstates at 
the detectors, and can only measure averaged oscillation. 
This loss of coherence by the detector during measure- 
ment is equivalent to incoherent v mass eigenstates ar- 
riving at the surface of the Earth (see, e.g. Ref. 13|). 
The resulting total v flavor conversion probability at a 
detector is 



Epsrc pffi 
i 



(5) 



Here P src and P® are conversion probabilities from the 
flavor to mass and from the mass to flavor states, re- 
spectively. The PMNS mixing matrix for the flavors and 
mass eigenstates are denoted with U, and A® is the tran- 
sition amplitude after propagation inside the Earth. The 
conversion probability at the source is neglected in the 
second equality, which is valid because of a negligible 
matter potential at the v production site. (However an 
interesting situation may arise if i/'s pass thr oug h the 
RG, which will modify the source probablity.) [23[ 

For the position of V407 Cygni, Right Ascension 
(RA) = 315.55° and Declination (DEC) = 45.74°, z/s 
travel a path length L — 2R^ cos8 n inside the Earth, 
with radius R®, and a nadir angle n = 90° — DEC. 
We calculate A® with a numerical code |l4T| that uses 
the Preliminary Reference Earth Model [15j for the den- 
sity profile inside the Earth, and v mixing parame- 
ters: Am\ x = 2.4 • 1(T 3 eV 2 , Am\ x = 8 • 1CT 5 eV 2 , 
sin 2 6x2 = 0.31, 6» 23 = tt/4, sin 2 13 = 0.02 and the CP 
violating phase <5 = 0. We consider normal v mass hier- 
archy only. The total conversion probabilities are plot- 
ted in Fig. [5] for neutrinos. The upper and lower panels 
correspond to v e — > v a and — > v a conversions, respec- 
tively. The anti-neutrino conversion probablities are not 
affected by matter in case of normal v mass hierarchy 
and vacuum conversion formalism apply. 

The most prominent features in the probability curves, 
the dip in the P„ c _s.„ c (and associated peaks/dips for 
other i/'s) at ~ 6 GeV, can be understood analytically 
from 2-flavor v oscillation framework with 1-3 mixing 
16]. For 9„ > 34°, which is the case for V407 Cygni, i/'s 
do not pass through the Earth's core. Conversions mostly 
take place in the mantle with an average density of (p) ~ 
5 g cm -3 . For this density, the low MSW resonance en- 
ergy is E L = Am 2 12 cos26»i 2 /[2v / 2G F (p)] « 0.08 GeV and 
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FIG. 2: Flavor conversion probabilities [Eq. for neutri- 
nos created at V407 Cygni and detected by an Earth-based 
detector at the South pole. The v £ — > v a and — ¥ v a proba- 
bilities are plotted in the upper and lower panel, respectively. 
We used normal hierarchy of v masses and a vanishing CP 
phase along with the best-fit oscillation paramaters (see text 
for details). 



1-2 mixing is strongly suppressed for our energy range of 
interest, E v > 1 GeV. The dip at ~ 6 GeV for P„ e -^ Ve m 
Fig. [2] (upper panel) corresponds to the 1-3 or high MSW 
resonance energy Eh = Am 2 3 cos 28 13 /[2\/2Gf(p)] ~ 6 
GeV. The width of the dip is 2t&n29 13 E H ps 3.6 GeV. 
At energies 3> Eh , the conversion probablities are dom- 
inated by vacuum oscillation. 

The v fluxes at a detector buried under the Earth's 
surface are calculated from Eqs. (HI and ([S]) as 



$?, rc p, 



(6) 



Event rates in a v detector depend on these fluxes and 
we calculate the expected number of events from a 7-ray 
nova outburst such as Nova 2010 in V407 Cygni next. 



IV. HIGH ENERGY NEUTRINO DETECTION 

The total number of neutrino-nucleon charge current 
(cc) interactions by u a , which produce secondary leptons 
of type a that may be detectable, can be calculated as 



W, 



N T t 



dE v J drt v eS (E v ,n) 



Vdet 

x [o?(E v )$%(E v , 0) + o?(E v )$*?(E v ,Slj\ , (7) 

below the energy (~ 100 TeV) at which v absorption in- 
side the Earth is negligible (see e.g. [Hj])- Here t is the 
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duration of the v outburst, Vdet is the instrumented de- 
tector volume and iVx is the total number of nucleons in 
that volume. The effective detector volume Kff ~ Vict 
for contained events at low energies. At very high en- 
ergies V c ff > Vdet, specially for because secondary 
/x's that are produced far from the physical detector can 
propagate inside due to their long range. We parameter- 
ize the cc cross-sections in the ~ 10 - 100 GeV range [ij} 
as 

< c = 7.3 • 10 _39 (J^/GeV) cm 2 
of = 3.8 • 10- 39 (-E„/GeV) cm 2 . (8) 

For a detector such as the IceCube Deep Core [Isj . 
with a 10 Mt fiducial volume, the total number of nucle- 
ons inside Vdet is Nt — 10 13 Na, where Na is the Avo- 
gadro's number. All astrophysical sources are essentially 
point sources in u's and one needs to carry out integration 
over E v only in Eq. (J7J. The lower range E v ^ depends 
on the energy of the secondary lepton at the detection 
threshold. The average muon energies are ~ 0.52i?„ for 
neutrinos and ~ 0.66-E„ for anti-neutrinos in the range 
E v ~ 10 - 100 GeV. For a rough estimate of N Va we as- 
sume a common E Utt h for both z/s and z?'s for the same 
lepton energy. Here we consider v^s and v^'s only, as 
the detection prospect of ~ 5 GeV muon tracks is better 
than the showers created by other flavors. For t = 15 
day, the duration of the 7 ray outburst from the nova in 
V407 Cygni, and for V c g = Vdct we calculate 

N V/t « 0.5±£J , w 0.3±g;l ; 10 < E v /GeV < 100, 

using the ir° model for 7 rays with the best-fit parame- 
ters, and their ±la variation. 



V. DISCUSSION AND IMPLICATIONS 

We have calculated ~ 0.1-100 GeV v fluxes from the 7 
ray Nova 2010 in the symbiotic binary V407 Cygni, from 
7T decays while assuming that the 7 rays detected by 
the Fermi LAT are from 7r° decays. Detection of these 
v's in the ~ 10-100 GeV range depends critically on the 



detector's angular resolution and uncertainty in the at- 
mospheric v background. For the IceCube Deep Core 
the atmospheric + events during 15 days of the 
7-ray nova is - 60 for AO = 10° following Ref. [3]. 
Future large scale detectors with a better angular reso- 
lution, such as ~ 3° in the Super-Kamiokande fl9j |. can 
reduce the background roughly by an order of magnitude 
to improve the signal-to-noise ratio, as well as measure 
background more precisely to reduce uncertainty. More- 
over, a similar 7-ray nova in the nearest symbiotic nova 
system AG Peg (D ~ 0.6 kpc), among a dozen known, 
would improve the signal by a factor of ~ 20. More en- 
ergetic 7 ray nova is also possible. The kinetic energy 
release in 2006 nova outburst in a nearby (D ~ 1.4 kpc) 
symbiotic system RS Ophiuchi is inferred to be an or- 
der of magnitude higher [20( than in V407 Cygni 2010 
nova 

The rate of classical novae in our galaxy is ~ 30 year -1 
in all binary systems (2lj . The rate of novae among ~ 
100 known symbiotic binary systems (22J, which seem to 
require for particle acceleration and 7 ray production, 
is more uncertain, once every few years. Following up 
of future novae with v detectors may establish a rate 
independent of the 7 ray instruments and answer the 
question of a hadronic or leptonic origin of the 7 rays from 
them. Non-detection of z/s from 7 ray novae can be used 
to constrain the emission models and system parameters 
as well. 

While astrophysical 7 ray sources have been hypoth- 
esized as potential TeV - PeV v sources, mostly from 
it's and if's created in pj interactions, relatively low- 
energy (~ 10 GeV) v sources which are simultaneous 7 
ray sources were thought to be rare. Surprising discovery 
of the nova in V407 Cygni by the Fermi LAT adds a new 
v source candidate class where, like supernova remnants, 
pp interactions may produce both 7 rays and z/s. 

We thank C.C. Cheung, CD. Dermer, A.Yu. Smirnov, 
C. Spiering and K.S. Wood for helpful comments and 
discussion. Work of SR is supported by NASA Fermi 
Cycle II Guest Investigator Program and the National 
Research Council. 
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